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Abstract 

I We consider the Gibbs representation over space-time of non-equihbrium dy- 

' namics of Hamihonian systems defined on a lattice with local interactions. We 

. first write the corresponding action functional as a sum of local terms, defin- 

ing a local action functional. We replace the local system by a translation- 
' invariant system whose dynamics has an identical space-time characterization. 

We study in details the irreversible properties of the new dynamics, define the 
local conductivity and show its equivalence with the Green-Kubo formula. 
Given the definition of the local heat conductivity and using conservation 
cp ■ of energy, we derive the shape of the temperature profile. Next, we find an 

explicit formula for the non-equilibrium stationary measure of harmonic sys- 
Tjj- I tems. Finally, we apply our scheme to various approximations of anharmonic 

■ Hamiltonian models, show how to compute their thermal conductivity and 

. recover results confirmed in numerical simulations. 
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1 Introduction. 



Non-equilibrium stationary states (NESS) of systems of classical Hamiltonian os- 
cillators located on a lattice coupled to heat baths at their boundaries have been 
extensively studied over recent years [71[9l[2l]. Numerous global results concerning 
those NESS have been obtained, including existence and uniqueness of the station- 
ary probability measure. The positivity of the entropy production and the validity 
of the Gallavotti-Cohen fluctuation theorem has also been established P, [2DI ES]- 
They have also been abundantly studied with the help of numerical simulations, see 
[5], [T7] for a review of those and a general overview of the subject. 

However, in contrast to equilibrium states for which the explicit Gibbs formula 
may be used, the effective computation of correlation functions of the stationary 
states out of equilibrium remains a challenging problem. Conceptually, one may 
distinguish two different reasons for that. The first one is common with the equilib- 
rium situation and has to do with the difficulty of dealing with nonlinear interactions 
between the components of the system. The second one has to do with the lack of 
translation invariance which is, in a sense built-in in the non-equilibrium set-up. 
This paper adresses the second issue for Hamiltonian systems defined on a lattice 
with local interactions. The difficulties we have just mentioned are particularly obvi- 
ous when one tries to identify the physical mechanism giving rise to a finite thermal 
conductivity in non-linear Hamiltonian systems and compute its dependance on the 
microscopic interactions. The conductivity is a local property of the system that 
relates local quantities like the current and the local temperature gradient. The goal 
of this paper is to outline a method which, by dealing with the lack of translation 
invariance of the system, allow to focus on the effect of nonlinearities on the non- 
equilibrium steady states and perform concrete computations of their correlation 
functions. The method is a priori independent of any approximation scheme and 
apply to any local Hamiltonian lattice system out of equilibrium. 

As we will argue later, what matters locally and away from the boundaries in a 
non-equilibrium set-up is not so much the spatial geometry but rather the specific 
way in which the time reversal symmetry of the equilibrium dynamics is broken. 
It turns out that heuristically, the non-equilibrium dynamics may be characterized 
locally in a simple, universal and translation-invariant way, see fl25l) below. This 
follows from local (on the lattice) conservation of energy and the characterization 
of NESS in terms of probabilities over pathspace introduced by Maes pLSj. There- 
fore, in order to study the local dynamics, we will consider a spatially homogeneous 
Hamiltonian chain of oscillators coupled to a stochastic thermostat at a fixed tem- 
perature T and, in order to recover the time-reversal symmetry breaking induced 
locally by the heat baths located at the boundaries, we include an additional non- 
Hamiltonian term in the deterministic part of the dynamics. The construction of 
the dynamics is very similar in spirit to the Evans heat flow algorithm [10] but 
ours proceeds through the introduction of a stochastic thermostat, which makes 
analytical computation straightforward. In particular, the proof that the dynam- 
ics satisfy the Gallavotti-Cohen fluctuation theorem is simple and the equations 
for the correlation functions in the stationary state may be easily written down. 
Systems with Langevin stochastic thermostats and non-conservative forces produc- 
ing entropy have been considered in [T3l [T9] . Also, when comparing our dynamics 
with the Nose- Hoover type thermostats (i.e. purely deterministic), it is important 
to notice that the additional deterministic term is itself time-reversible. It is only 
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when coupled to the Langevin thermostat that the dynamics becomes irreversible in 
time, in a stochastic sense. We shall see how to recover the linear response theory's 
Green-Kubo formula in a natural way. However, our scheme hopefully also allows 
to describe the behaviour of systems far from equilibrium. We shall also see how 
to reconstruct the temperature profile on the original chain coupled to heat baths 
at different temperatures. The basic assumption that is made is that two systems 
having the same time-reversal symmetry breaking functional will have basically the 
same average heat current. The validity of this assumption would certainly benefit 
from more of theoretical and experimental arguments. 

In the last section, we consider a concrete anharmonic Hamiltonian model, 
generic of models having a normal heat conduction behaviour. The Hamiltonian 
is given by, 

H{q,p) = E [| + + \qt + Y(g. - <lr^if] (1) 

On that example, we will see how to implement our strategy and study the prob- 
lem with the help of the non-equilibrium translation invariant dynamics we have 
defined before. We will show that this dynamics provides a source term, absent in 
equilibrium, in the equation for the evolution of the average heat current j in the 
chain. This is of course what should be expected because the dynamics is built so as 
to break the time-reversal symmetry of the equilibrium dynamics in a proper way. 
Roughly speaking we get an equation of the form, 

j^m = -im+mt) + s{t). (2) 

The first term on the RHS of this equation is a damping term coming from the 
coupling of the chain to an external friction represented by the coefficient 7. The 
second term comes from the Hamiltonian evolution and vanishes when the anhar- 
monic part of the interaction is zero, i.e., when A = 0. The third term represents 
an external "creation" of current and comes from the action of the non-equilibrium 
dynamics. Equation ([2]) will therefore yield a balance equation for the heat current 
in the stationary state. First, we review the linear case (A = 0), for which we give 
an explicit formula for the stationary measure out of equilibrium. Next we study 
two approximations which have been used recently to analyze the effect of the an- 
harmonic terms on the non-equilibrium properties of the Hamiltonian system. The 
first one [2[ |12] amounts to replace the effects of the anharmonic terms by a ran- 
dom exchange of energy between nearest-neighbours, this gives rise to a normal heat 
conductivity and to a linear profile of temperature. 

The second one has been derived and used over the years since Peierls P, HH 
[HI [151 [22l [2^- It amounts to a closure assumption on the hierarchy of correlations 
and to a derivation of a stationary Boltzmann-type equation with a collision ker- 
nel representing the interaction between the phonons due to the nonlinearity of the 
system. The basic mechanism which gives rise to a normal heat conductivity has 
been known to physicists for a long time, it was postulated to be the result of some 
particular resonances in the collision between phonons. However, the lack of explicit 
localization and of a proof of the existence of those resonances prevented the effec- 
tive computation of the thermal conductivity in terms of the physical microscopic 
parameters. The localization of the resonances in the model ([T]) for large pinning 
has been achieved in |15J, where a similar translation invariant dynamics has been 



devised, and the value 
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for the thermal conductivity was obtained. The value of a = 0.275637 was computed 
in [I] where the expression was also tested and numerically confirmed. The 
localization of the resonances at = was obtained by Pereverzev ^23j who used it 
to discuss anomalous heat conduction in a related anharmonic model. 



2 Local time-reversal in models for heat conduc- 
tion. 

In order to explain the behaviour of the thermal conductivity in crystalline solids, 
one often modelizes the solid by a chain (or lattice in higher dimension) whose ends 
are coupled to heat baths maintained at different temperatures. The coupling can be 
taken stochastic and more precisely of Langevin type. In one dimension, the set-up 
is as follows. At each site i of a lattice {1, . . . , N} is attached a particle of momentum 
Pi and position gj. The dynamics is hamiltonian in the bulk and stochastic through 
the Langevin coupling to heat baths at the boundaries. The Hamiltonian is of the 
form, 

N N 

H{P,q) = E(|P' + V{q^)) + E U{q^ - q^-l) + U{q,) + f/(g^). (4) 

i=l 1=2 

The equations of motions are given by, 

dqi = pidt, i = l,...,N, 
OH 

dpi = - — {p,q)dt, i = 2,...,N -1, (5) 



dqi 



dpi = -t:—{p, q)dt - -fpidt + \/2'y kTi dwi 
dqi - - V 

OH 



and. 



dpN = — (p, q)dt - 7PArdt + J 2'ykTii dwr . (6) 
oqN 

Tl and Tr stand for the temperature of the left and right reservoirs, respectively, 
whereas wi and Wr are two independent standard Wiener processes. 

It is an easy fact to check that when Tl = = T = (3~^, the equilibrium 
measure on the configuration space R^^ whose density with respect to the Lebesgue 
measure is given by 

p(p,g) = Z-ie-'3^(M) (7) 

is invariant (stationary) for the stochastic dynamics defined above. In the case of 
two different temperatures, existence, uniqueness and exponential convergence to an 
unique invariant state has been established under fairly general conditions on the 
potentials U and [TJ [9l [23] . In the case of harmonic coupling, the covariance of 
the stationary state has been exactly computed in [211 [25] . 

In order to describe the conduction of heat in the crystal, one defines a local 
energy function, 

2 

h^{p, q) = ^ + Viq,) + i(f/(g.+i - g.) + f/(g. - g.-i)) (8) 



for i^l.N, 

2 

hiip,q) = y + Viq^) + U{q,) + \U{q2 - Qi) (9) 

and similarily for h]\[- The heat current is defined through the time evolution of the 
local energy, 

(10) 

for i ^ 1,N and where, 

Jt = - (li+i){Pi+Pi+i), (11) 

ii = ^F{qi_i - qi){pi+Pi-i), (12) 

where F = —U' and ji is defined to be the microscopic current of energy or heat 
between atom i and i + 1, i.e the rate of energy transfer between those atoms. We 
observe that = j~ and in the following we shall use the shorthand notation, 

ji — Ji ■ 

We use the notation (.) to denote average with respect to the stationary state. And 
define the local temperature in the chain to be 

(P?) = (13) 

and the average heat current is, 

j = {ji) = Y ^^^^^ " Q^+l)iP^ + Pi+l)) (14) 

which must be constant throughout the chain due to conservation of energy. 
This follows from, 

=or)-a-)=o (16) 



dt 



and j + 1 = ji . 
Fourier law states that, 



J = /c(T,)(T,+i-T,) (16) 

where k is the conductivity of the crystal, a major problem is to understand how 
the anharmonic part of the Hamiltonian makes this constant finite and to compute 
its value from the microscopic parameters. 

An elegant and systematic way of describing non-equilibrium dynamics was de- 
vised by Maes [18]. For stochastic processes, it amounts to compare the probability 
of trajectories of the non-equilibrium process at hand and of its time-reverse. By 
the Girsanov formula, the ratio of the two probabilities is given by the exponen- 
tial of an "action functional". Then the part of the action functional which breaks 
the time-reversal invariance is identified with the entropy production. Let the 
pathspace measure on trajectories ^ = {{q{s),p{s)), s G [— t,t]) obtained from the 
dynamics (15|6|) started from initial conditions which are distributed according to 
some measure /i. We define the time-reversal operator on trajectories as follows, 
(n^)s = (g(— s), — p(— s)) and let P^II be the pathspace distribution of the time- 
reverse process. In [201 126] , it has been show that the ratio of the two measures may 
be written as, 

r/P* 

^ = exp(P*(0) (17) 



modulo temporal boundary terms involving the distribution /x. 
The action functional i?* is given by, 



+ / -jN-l{s)ds - (/lAr(^t) - hN{^~t)), 



with = Tj^^ and Pr = T^^. We observe now that for any collection Ki, i 
1, . . . ,N, such that Ki = and K]^ = Tji, 



N-l t N . 

1=1 * i=l * 

where we defined the local entropy production cTj, 

= - -f^)J^■ (20) 

This follows from local energy conservation inside the chain written under the inte- 
gral form, 

hki^t) - hk{i-t) = j\jk{s) - Jk-i{s))ds (21) 

Except for the arbitrariness of the -ft^j's, this shows that the conservation energy 
imposes strict conditions on the representation of the action (fT8|) in terms of the 
local variables ji and hi. 

Among all the equivalent representations for the action difference ( lT9l) . we now 
choose one which is useful because it allows detailed local analysis of the dynamics. 
Let us take Ki = T{^), with T(.) a smooth function over the interval [0, 1]. Now, 
at lowest order in A^~^, ffTOl) becomes. 



1 ft \/T(^) ^ ?■ 
^*(0 = T7E/ Ms)^f;^ds-j:Pi^)mt)-k{^^,)). (22) 

In the limit iV — >■ oo and around a fixed point x of the interval [0, 1], the variation 
of the coefficient in front of the integral of the current and the local energy goes 
to zero. So we subdivide the interval [0, 1] in subintervals of size e, with e small 
but N >> e~^, each interval being centred about some point that we denote by Xk- 
Now, at lowest order in e, we write, 

R\0 = T.RW (23) 

k 

The local action functional -R^(0 is given by, 

^^(0 = WT^ E f Ms)ds-(3{x) E (h^i^t) - h,{^.,)) (24) 

where B^{x) = {j\ \j — Nx\ < ^Ne}. Note that |-Be(x)| = eA^. So locally around the 
point x, we get the following characterization of the non-equilibrium dynamics. 



^- = ^^pfw^^ ^ f Ms)ds-(3{x) E mt)-hM.,))] 



(25) 



with, 

^=n«« (26) 

The important feature of (!24|) is that the original problem has become locally (in 
the box B^{x)) translation invariant. For a given x, we now fix T = T{x) and 
r = N^^VTlx), M = [Ne] and in the next section, we define a dynamics which 
has exactly the property (1251) for those parameters, see (jSHD- Then, in the last 
section, we will write the equations that the correlations of the stationary state 
satisfy in a concrete model. This amounts to ask local stationarity in the original 
non-translation invariant chain and allows on one hand to interpret the parameter 
T as a local temperature because it gives the average kinetic energy. On the other 
hand, it fixes the current as a function of the temperature gradient VT and local 
temperature T. Going back to the original problem, the fact that the average current 
in the stationary state is constant in position then fixes the kinetic energy (i.e. 
temperature) profile as we shall see at the end of section 5. We have derived (!25l) for 
Hamiltonian containing only nearest-neighbour interactions. Presumably the same 
identity may be derived for systems with sufficiently local interactions. And in the 
next section, we devise a dynamics which include more general interactions. 



3 The dynamics and its entropy production. 

We consider a periodic lattice hamiltonian system described by the Hamiltonian, 



M 



H{q,E)=Y: 



1=1 



Pi 
2 



M-l 



fc=i 



(27) 



with the convention i + M = i. Note that the form of the Hamiltonian is slightly 
more general than in the previous section but is the same when U'^ = 0, for k ^ 1. 
Proceeding exactly as in (|8|T0l) . one computes the local energy variation rate, and 
by analogy, 

Ji = k T.iPi+Pi+k)F''{qi - qi+k), (28) 



2 

k 

where = —{U^)' . ji is interpreted as the total heat current entering the site i 
coming from the "right" of the chain. We shall use below the spatial average of the 
energy current, 

^ M 

J=jfEl. (29) 

Our dynamics is, 
dqi = pidt 

dpi = --fPidt - —dt + ^(F^(gi_fc - g^) + F^{qi - qi+k))dt + ^2^dwi 

(30) 

and the Wi are standard independent Brownian motion i = 1, . . . , M. The term 
proportional to r is the non-equilibrium part of the dynamics and is responsible for 
the breaking of the time- reversal symmetry of the equilibrium dynamics (at r = 0). 



7 



We will see that its particular form allows to show the validity of fl25|) and thus the 
validity of the fluctuation theorem and of the Green-Kubo formula for the thermal 
conductivity. The generator of the dynamics is, 

L = Lo + Lr (31) 

with, 

L^^y _^p^A_ _ ^ A + ^ A + ^T— (32) 

° i 'dpi dqi dpi dpi dqi dp] 

and 

^r = ^ E(^'(?-^ - + - (33) 

i,k 



We note the basic identity, 



L^rH = UH = ^ (34) 



where is the adjoint of the operator L^. with respect to the Lebesgue measure. 
We assume the existence, uniqueness, smoothness and regularity in the parameter 
r of a stationary measure satisfying, = 0. 

Remark. The structure of the additional force proportional to the parameter r is 
chosen so that is satisfied, or more generally so that the generahzed detailed 
balance relation holds, see (152|) below. This in turn implies the validity of the time- 
reversal characterization of the non-equilibrium dynamics (l25l) (and ( |56i) . see below). 
The only degree of freedom left is the value of the friction parameter 7 which does not 
enter in fl25l) and fl56|) as it does not in the stationary measure at equilibrium either. 
It should be seen as a regularizing parameter. Indeed as we shall see explicitely in 
the last section, its main role is to provide an "external" degradation mechanism 
for the energy current in the chain. The physical degradation should come however 
from the nonlinearity of the chain and we shall take the friction 7 — in order 
to evaluate the conduction properties of the chain. But before taking the limit 
7 — > 0, we will always take the limit M ^ 00. By doing this, we will be able 
to discriminate between normal heat conduction and anomalous conduction, in one 
case the conductivity is infinite in the other it is not. It is in the limit of 7 — that 
we should expect the system defined by the dynamics (15(1 and the local dynamics 
of the original non-equilibrium problem to be equivalent. 

We first study the evolution of the Shannon entropy under the evolution defined by 
the generator L. To any measure with a smooth density / on the phase space R^^^ 
we associate the Shannon entropy, 

Sif) = - Jdxfhgf (35) 

with the notation x = {p,q)- Compute its derivative with respect to time for a 
density whose evolution is given by 

dtft = L^ft (36) 

Next, 

dtSift) = -j dxil + log ftM ft) (37) 
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J dxil + log ft)L'' ft (38) 
dxftLlogft (39) 



dtSift) = -J dxftL{\og ^)+pJ dxftLH (40) 



ft 

where p° = Z^^ exp{~-(3H). For the first term, we note the identity, 

L(log h) = h-^Lh - h-'^r{h, h) (41) 
where F is the Dirichlet form associated to the generator L, 

So that (HOj) becomes, with the notation gt = ft{p^)^^ 

dtSUt) = - j dxp'^Liip'^r'ft) + J dxp'^g^'Tigt, gt) + P J dxftLH (43) 
Integration by parts of the first term yields, 

dtSift) = -j dxiip')-'ft)L^p' + I dxp'gi^T{gt, gt) + P J dxftLH (44) 
Using flMl) . one gets L^p" = p'^H^ and therefore, 

dtS{ft) = -j dxft^ + J dxp'gt^T{gt, gt) + [3 J dxftLH (45) 

When ft is the density of the stationary measure p"^, the last term drops and one 
gets. 



= - y dxp^ix)^ + J rfxp"(x)F((p")-V^ (46) 
It is natural to define the entropy production as 

^ = ^ (47) 

And using the fact that the second term of fj46l) is given by a (positive) Dirichlet 
form which vanishes iff p^ = p°, one concludes that the average entropy production 
in the stationary state is positive and that the current J 7^ iff r 7^ 0. 



4 Generalized detailed balance relation. 

4.1 Fluctuation theorem for the entropy production. 

Consider the ergodic average of the entropy production, 

1 If* If* Mr 

at ^ -W{t) ^ - a{s)ds = - y —J{q{s),p{s))ds. (45 



If the process is ergodic one has, 



hm at = {a) (49) 

t — ^oo " 

for p^-almost every initial conditions g(0),p(0). The fluctuation theorem is con- 
cerned with the specific shape of the large deviation functional e{w) of at defined 
roughly as, 

P(o-t e [w - e, w + e]) ~ exp(-e(w)t) (50) 
The Gallavotti- Cohen fluctuation theorem states that 



e{w) — e{—w) = —w. (51) 

and by standard arguments [H [131 [HI [26] , assuming good properties of the process, 
the validity of this identity is basically ensured by the generalized detailed balance 
relation which is elementary to check in our case. Indeed, defining the reversal of 
velocities operator vr, 7r/(p, g) = f{—p,q), one checks, for any a, 

(p°)-V(L^ - aa)7rp° = L - (1 - a)a. (52) 



4.2 Time-reversal and entropy production. 

We now show that the dynamics f[30l) has the property f[25|) . expressing the fact the 
time-reversal symmetry is broken in a proper way. In order to compare the space- 
time probability of a process and of its time-reverse, we first choose a reference 
process, to which we shall compare both. Let the pathspace measure on trajec- 
tories ^ = ((g(s),p(s)), s G [— t,t]) obtained from the dynamics ( l30i) started from 
initial conditions which are distributed according to some measure p. We define the 
time-reversal operator on trajectories as follows, (n^)t = 7r^_t, where vr as above 
reverse the sign of the momenta. We want to evaluate the ratio dPji/dPlJl, where, 
as in section 2, P^II is the distribution of the time-reverse process. We compute 
first dPji/dPpO, where Pjo is the pathspace measure obtained from the dynamics 
( 1301) with r = and as initial distribution the Gibbs measure p° = exp{—j3H). 
A direct application of Girsanov formula yields. 



dPl = exp{A\0)K'> (53) 



where, 



+ £ -^F{q{s))p{s)ds + Inp(e-t) - \np\^^t), (54) 

with, Fi{q) = J2k{F''{qi — qi-k) — F^{qi — qi+k)) and the first term of (1M|) is an Ito 
integral. Using the reversibility of the equilibrium measure, i.e. P^o = -PpoII, we 
write. 

Being invariant under the time reversal operator 11, the second term in ([M|) does not 
contribute to the difference in the above f[55l). The fact that the first term does not 
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contribute is somewhat less obvious, see for instance example 3 in ^19j for a proof 
of this fact in a similar set-up. Finally, using F{q{s))p{s) = 2MJ{s), we get, 

R\0 = A\0 - = j\ <s)ds + mi-t) - mit). (56) 

modulo the boundary terms involving the distribution /i, this is the same as (125!) 
with the identification r = iV-^VT(x), jS = T-^{x) and M = [eN]. 



5 Properties of the stationary state. 

5.1 Heat conductivity and Green-Kubo formula. 

We define the heat conductivity k as 

K = lim lim limr^^ (J),t (57) 

and show now that it formally coincides with the usual Green-Kubo formula for 
Hamiltonian systems. The Green-Kubo formula for the thermal conductivity of an 
Hamiltonian system is, 

M 



KGK = lim — / {J{Q)J{s)),ds. {hi 



where J(s) is given by (1291) in terms of the coordinates q{s),p{s), solution of the 
Hamiltonian deterministic equations (i.e fl30|) with r = and 7 = 0) with initial 
conditions g(0),p(0). 

As before, we assume the existence and uniqueness of a smooth (both in space and 
as function of r) invariant measure for the process fl5U]) . whose density p"^ satisfies, 

LV = (59) 

Developping at first order in r the density = + rp^ + . . ., we get an equation 
for the first-order correction to equilibrium (at temperature T) from, 

L^p^ = L^p° + tL^p^ + LIp"" + . . . (60) 

and thus, formally, 

P' = -{Llr'Llp' (61) 
We first compute limT-^o {JY ■ Using (I6T1) and again (IMI) 

= -{LIVy-Jp'. (62) 



Next, because the expected value of the current is zero in equilibrium and ((Lq 



T\-l\T 

sc Liic c^pcuicu vaiuc ui iiic uuiiciiL lb Ziciu Hi cquiiiuiiuiii aiiu yyu 

Lq \ we get 

I TV M /V /"OO rr. 

lim^ = -^(J(Lo)-V)^o = {JmJ{s)y'^)^ods (63) 

where (J(s))''^'^ is the solution of 

-f(J(.)P^ = Lo(J(.)r^ (64) 
ds 
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with initial conditions (J(O))'^' = J(0). Formally, the most natural thing to do is to 
take next the limit of M — > oo and then the limit 7 — to recover the Green-Kubo 
formula (I58p . So that we get finally our definition of the conductivity ([57]). 

Remark. The discussion of these limits is a bit delicate in the case of anomalous 
conductivity in which case one tries to evaluate the dependance in M of the con- 
ductivity K, proceeding as we do yield k, = 00 each time that there is anomalous 
conductivity. If one wishes to study the size dependance of the conductivity as for 
instance in t3ii23j, it seems that one should take 7 scaling as and let go M — cx). 
But at the moment, we do not have any justification for that particular scaling of 
the friction and suspect that it might be model-dependent while the definition (1571) 
is not and does discriminate between normal and anomalous heat conduction. 

5.2 Temperature profile. 

We come back now to the original problem of the chain heated at its boundaries. 
As we have seen above, if the conductivity is finite, one may write, 

{jy = k{T)t (65) 

where k will in general depend not only on the temperature T but also on the 
physical parameters of interaction, in the limit of large N and small 7. As we shall 
illustrate in the next section on a concrete example, one obtains such a relation by 
simply writing down the stationary Fokker-Planck equation for the dynamics (1301) 
and analyzing the effect of the nonlinearity on the nonequilibrium properties. But 
now, as we explained at the end of section 2, we do the identification r = N~^VT{x), 
T = T(x), M = [eN] so that ([56]) is really identical to ([25]) . And we get, 

{J^) = {jy = KiTix))^VTix). (66) 

for i G B^{x), the box of size e around the point x. This is where one explicitely 
uses the assumption of equivalence between the original local dynamics (in the box 
B^{x)) and of the dynamics ([3D]) in the limit 7^0 because of their identical time- 
reversal property (!25l) . ([56]) . As we observed in ([T5]) . in the stationary state, the 
energy current must be constant along the chain and therefore, one gets, 

Vk{T{x))VT{x) = 0, (67) 

away from the boundaries. The temperature profile is finally fixed by imposing the 
boundary conditions, which also sets the value of the current as a function of the 
imposed temperature difference. 

6 Models for heat conduction and current degra- 
dation. 

Having investigated the general properties of the stationary measure of the dynam- 
ics ( 1301) and shown that its spacetime extension has the natural characterization 
(1251) . (ISBI) . we now study the stationary measure of a concrete model by writing 
down the equations that the correlation functions satisfy. As a generic model for 
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anharmonic systems having a normal conductivity, we consider a lattice system 
described by the Hamiltonian, 

M 2 ^2 \ , ,2 

H{%P) = E [| + + -At + Y^Q^ - Q^-if] (68) 

i=l 

and the non-equilibrium dynamics fl30l) with periodic b.c. becomes, 
(igi = pidt 

dPi = -IPidt - a;^((2 + - - g^+i) - Agj' 



+ ^(g^+i - + plTdw, (69) 
As above, one defines a local energy function, 

K{p_.q) = I + + + ^(g. - g.-i)' + {q^ - q.+if (70) 

and the heat current in the chain ji, 

Ji = Y(?i+i - Qi){Pi+Pi+i)- (71) 
and because of translation invariance, we shall often use the notation, 

J ^ iW ■ (72) 

for the average microscopic current in the stationary state. Compute now the evo- 
lution of the current ji under the Hamiltonian evolution, 

= y((p?+i-pD-(?»+i-?i)^^((2+/i)(gi+gi+i)-(gi+i+g*+2)-(g*-i+g*))+A(g^'+gf_^ 

(73) 

As the dynamics is translation invariant, we can assume that the unique stationary 
measure is also translation invariant, which implies that, 

\ dt , 

One obtains in the stationary state, 
d ... /d^jiV 



^{{q^qhl) - (74) 



(J.) = ( ^ ) - 7 {J^r + ^ ((g.+i - q^^.ry = (75) 



dt ' \ dt I ' ' 4T 
So that we finally get 



4 

- 7 {3^Y + ^ ((g.+i - q^-l?Y + A {{q^qUY - (qh+iY) = (76) 

This is the basic equation for the current. The first term is a damping term for 
the current, the second term is a source term for the current (as we see below) and 
has the same sign as the parameter r, it comes directly from the irreversible part 
of the non-equilibrium dynamics, i.e it is produced by the force proportionnal to r 
in (1301) . The last term comes from the Hamiltonian part of the dynamics. It should 



be responsible for the "internal" degradation of the current in case of a normal 
conductivity. 



Harmonic case. 

We look first at the harmonic case (A = 0). By a direct computation (i.e. by 
checking L^p = where is the adjoint of the operator defined in ( !3T]) ). one can 
show that the measure over the phase space R^*^ with Gaussian density, 

p{q, p) = exp (^-(3H{q, p) + ^J{q, p)j (77) 

is stationary for the process defined by fl69l) when A = 0. We now look directly at 
the equation fl76l) with A = and at the lowest-order in r since we are interested in 
computing the heat conductivity, i.e we write. 



7(j.r + ^((gm-*-if) =0 (7J 



4T 

where denotes the expectation value with respect to the equilibrium harmonic 
measure. This last term is readily computed and (ITHI) yield, 

uP'T r\ sin^(27rx) , , 

^^^y = / . 2^/- 2 . dx (79) 

7 J-\ jj"^ + A sm [TTX) 

for large A^. Here the basic problem of harmonic chains is apparent, when the exter- 
nal friction 7 is taken to zero the current becomes infinite, the harmonic chain has 
no internal mechanism which could degrade the current that is continuously created 
by the "thermal force" proportional to r (i.e by the local temperature gradient). 
According to our definition, the conductivity of the harmonic model is, 

K = lim lim lim r^^ (JV = +00 (80) 

Let us see now how to deal with the anharmonic interactions. Observe first that 
if the stationary measure was Gaussian, the conductivity of the system would be 
infinite, indeed under Gaussian assumption and translation invariance, 

{{QiQliY - {qh^+ly) = K{m^+lY (q-Y - {q^q^+lV (qliY) = o (8i) 



and one gets (ITHI) again. Observe also that (^qiqf^i^ — {q^qi+iY is odd under the 
exchange of the indices i and i + 1, as is the current, by definition ( |7T1) . 

Stochastic approximation of the non-linearity. 

The simplest possible assumption which can be made about the stationary state is 
to assume, 

{q^(i^.+lY - {qfqi+iY = (82) 

for some rj. The important point is that such a term may be produced by a dynamics 
which amounts to randomly exchange energy between nearest-neighbours in the 
chain and not by an external damping term in the dynamics. This random exchange 
is responsible for the degradation of the current constantly produced by the "thermal 
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force" proportional to r. See [2] for details and rigorous results on the model. Now, 
(176|) become 



- 7 {J^y - V {J^Y + ^ - = (83) 



AT 

and the conductivity is, 

oj"^ rk sin^(27rx 



-dx, (84) 



T] 7-i /i^ + 4 sin (ttx 

when /i = 0, one gets = ^ which is the number obtained in pj^ . As the conductiv- 
ity does not depend on the temperature T, from fl67j) . one gets a linear profile. More 
sophisticated models and stochastic approximations have been devised to tackle the 
problem of anomalous heat conduction [3]. 

Closure approximation. 

We refer to [15] for details about the closure approximation on the model we are 
considering in this section, but we give the outline of the argument. Observe that 
(1761) involves a four-point correlation function. In order to interpret it as term 
providing a internal mechanism for the degradation of the current, one would like 
to show that it is simply proportional to the average current, namely a two-point 
correlation function. To achieve this, the simplest thing to do is to assume that the 
stationary measure is Gaussian. As we observed above, this obviously factorizes the 
four-point correlation in terms of the two-point correlation but does not produce 
any significant correction to fl7^ with respect to the harmonic case. The next 
simplest thing to do is to write down the hierarchy of equations satisfied by the 
n-point correlation functions and assume that the six-point correlations factorize 
in terms of the two-point correlations. After linearizing around the equilibrium 
harmonic solution, this yields a linear expression of the four-point correlations in 
terms of the two-point correlations. The main point of the analysis is to identify 
the resonances in the interactions between the harmonic degrees of freedom (the 
phonons). To formulate the end result of the closure approximation we need to 
consider the general set of correlation functions, 

Ji = {piqi+iY (85) 

the fact the RHS does not depend on i follows from translation-invariance. In 
particular, Ji = j. We denote by J the vector made of the J^'s. Now, after closure, 
the stationary equation for this vector takes the form, for large /i, 

-7l-^jCU) + (L = (86) 

UJ' 

Where £ is a linear operator which does not depend on any physical parameter 
entering the equations of motion. The components of the vector a coming from the 
thermal force are given by, 

n /"i sin(27rx) sin(27rmx) , 
J -I fi^ + 4:Sm{TTX) 



^There is a slight difference coming from the fact that in [2J, the definition of the current 
contains a part coming from the random exchange itself, for small rj this becomes negligible in the 
computation of the conductivity. 



When computing a we have neglected correction of order A in the equihbrium mea- 
sure has been neglected. When 7^0, one inverts £ and gets, 

i=^^"^(-) (88) 
from which one may obtain Ji = j and thus the conductivity which is (for large /i) 

- = «^ (89) 

The value of a = 0.275637 was computed in where the expression ([3]) was also 
tested and numerically confirmed. As for the temperature profile, one gets, 

V^VT(x) = (90) 

away from the boundaries. The solution of (!90|) and the problem of its boundary 
conditions has been discussed in p. 

Finally, we note once again that our method is not limited to the example we 
are treating here. Each model will have a particular source term in the equation 
governing the degradation of the current. This source term can be directly computed 
from the dynamics (1301) . In particular, it would be very interesting to study the heat 
conductivity of the Toda lattice with this method, which we plan to do in a future 
publication. 

Acknowledgments. It is a pleasure to thank Cedric Bernardin, Thierry Bodineau, 
Alain Schenkel and Herbert Spohn for useful discussions. 



7 Appendix: Solution of the Harmonic case. 



We take the process defined by fl69|) with A = and solve the equations for the 
correlations in the stationary state. 

It is most convenient to work in coordinates where the equations become diagonal, 
and we introduce the Fourier coordinates for the periodic harmonic chain by 

1 ^ 

Q. = 47Ee^^'^'?. (91) 



with —N/2 + 1 < k < N/2. The Pk coordinates are defined in a similar fashion. We 
note for further purposes that = Q-k and Pk = P-k- In those coordinates, the 
Hamiltonian reads, 

N/2 |p2| 1^ 12 

H{Q,P)= E [^ + ^^^] (92) 

k=-N/2+l ^ ^ 

where 

ul = u\f^^ + 4sin\^)) (93) 
The spatial average of the energy current in the periodic chain fl29|) reads now 

1 2'Kk 

J=^ E sin(--)Im(Q_,Pfc). (94) 

k=-N/2+l 



In Fourier coordinates, the equations fl69|) with A = become, 

dQk = Pkdt (95) 
dPk = -iPkdt - {ool - I- sm{^))Qkdt + ,/2irdWk (96) 

The Wk are complex Wiener processes which satisfy Wk = W-k and, (Wk, Wi) = 
6{l + k). 

The equations for the correlations in the stationnary state are, 

{PkQ^k) + {QkP-k) = (97) 
Pk\') - ujl(m') = (98) 



27 (|P,|2) + tulak^ {PkQ^k) = 27T (99) 



27 {PkQ-k) - lu^iauT^ {\QkV) = (100) 

where = 2u;^ci;^^ sin(^). These are readily solved and yield 

(iPfcl') = 6kT (101) 

(PkQ-k) = h^'tSkakT (102) 



.,2 ^2 ^2. 



with, Sk = {I — -^^t) ^ ■ Keeping the lowest-order in r and doing a little compu- 
tation gives (1791) again. 
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